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ABSTRACT 

Pal 5 is a low mass, low velocity dispersion, globular cluster with spectacular tidal tails. We use the 
SDSS DR8 data to extend the density measurements of the trailing star stream to 23 degrees distance 
from the cluster, at which point the stream runs off the edge of the available sky coverage. The size 
and the number of gaps in the stream are measured using a filter which approximates the structure 
of the gaps found in stream simulations. We find 5 gaps that are at least 99% confidence detections 
with about a dozen gaps at 90% confidence. The statistical significance of a gap is estimated using 
bootstrap re-sampling of the control regions on either side of the stream. The density minimum closest 
to the cluster is likely the result of the epicyclic orbits of the tidal outflow and has been discounted. 
To create the number of 99% confidence gaps per unit length at the mean age of the stream requires 
a halo population of nearly a thousand dark matter sub-halos with peak circular velocities above 1 
kms -1 within 30kpc of the galactic center. These numbers are a factor of about 3 below cold stream 
simulation at this sub-halo mass or velocity, but given the uncertainties in both measurement and 
more realistic warm stream modeling, are in substantial agreement with the LCDM prediction. 
Subject headings: dark matter; Local Group; galaxies: dwarf 
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1. INTRODUCTION 



lAbelll (|1955l ) cataloged the globular cluster Pal 5 in the 
National Geographic Society - Palomar Observ atory Sky 
Survey, noting that Baade and iWilsonl (|1955[ ) had each 
independently discovere d the o bject , which Wilson listed 
as the Serpens cluster. lArd (119651) assi gned Pal 5 the 
lowest possible iShaplev fc Saw veil (|T?J2~ :■ concentration 
clas s. The cluster also has an unusually low luminos- 
ity (iSandage fc Hartwickl 19751). The quant itative King 
model (jvon Hoernerl 119571 : IKind IT9621 fl966h concentra- 
tion parameter is a notably, but not unique ly, low 0.52 
(jKing et al.111961 IWoltierlfl975l IHarrislll996l ). 

In the early photographic data the density profiles 
became uncertain toward the tidal radius as the clus- 
ter sinks into the distributio n of foreground and b ack- 
ground stars in the galaxy. iGrillmair et al.l (|1995l) ac- 
quired new, multi-degree, photographic data and intro- 
duced photometric selection procedures to screen out 
background stars, yielding the important discovery that 
near the tidal radius the stars were distributed in asym- 
metric clouds with a "striking resemblance" to the ex- 
pected tidal tails. The depth and ph otometric unifor- 
mity of the Sloan Digital Sky Survey (jYork et al.ll2000f ) 
have provided increasingly good views of the tidal tails 
of Pal 5. In the SDSS commissioning data the Pal 5 
tails were detecte d out to a distance of about ±1.3° 
from the cluster (Odcnkirc hen et ahl 120011) . Increased 
sky coverage of the SDSS and the introdu ction of a sta- 
tistic ally optimal matched filter technique (Rockosi et al. 
2002) allowed the trailing tail to be traced t o about 6.5° 
from the cluster ( Odcnkir chen et all 120031 ). Their re- 
sulting detection of density variations along the stream 
was a major new development. SDSS DR4 data im- 
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Fig. 1. — The match filtered star densities in the region of the 
Pal 5 stream in the SDSS A and 77 co-ordinate system. The raw 
image has been smoothed with a 3 pixel Gaussian. The object 
above the stream is the foreground cluster M5. 



proved both the sky coverage and the photometry, al- 
lowing [GrinmaiF^i^Ionatog ([200§) to trace the trailing 
tail to about 16° from the cluster and the leading tail to 
about 6°, where it runs off the edge of the current survey 
area. 

Current data shows Pal 5 to be one of the lower 
luminosity(Harris 1996) and lowest velocity disper- 
sion, l.l±0.2kms _1 (|Odenkirchen et al.l l20T)2l) . glob- 
ular clusters known. The l ong, thin, cool tidal tails 
(|Grillmair fc Dionatosl 120061: lOdenkirchen et al.l I2009D 
have significant density variations along their length. 
Consequently it has become a key system for the study 
of the small scale dynamics of the galactic halo with an 
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Fig. 2. — The matched filtered star map of the Pal 5 field, with Pal 5 and the foreground M5 cluster masked out. To remove the varying 
background, the masked image has been smoothed over 4° subtracted from the original image, and then smoothed with a 2 pixel, or, 0.2° 
Gaussian. The analysis is conducted on the original uncorrelated pixels. We have made no attempt to straighten the southern part of the 
stream, left of the cluster in this image. 



emphasis on the physics behind the origin of the star 

stream density variations. 

In a series of papers IKurroer et all lj2008_ 12010. l20ll 
use both dynamical analysis and precise n-body simula- 
tions to demonstrate that stars leaving a globular cluster 
at low velocity through the Lagrange points execute co- 
herent epicyclic orbits which lead to pileups at the low 
points of the cycloids and low density regions between. 
Their detailed n-body simulations show that the peri- 
odic structures are present for general cluster orbits, but 
that the range in phase angles emerging through the La- 
grange points cause the structures to blur out with dis- 
tance down the stream and that all quantities vary as 
the cluster and stream orbit around the galaxy. 

Dwarf galaxies also produce star streams although they 
do not have the two and three body interactions that 
provide a nearly continuous flow of stars to the zero en- 
ergy surface in a star cluster. The possibility of dynami- 
cal instabilities (|Comparetta fe G uillen 2011) that could 
create ga ps if they become nonlinea r has been raised. 
However, Schneide r fc Moorel (|2011D demonstrate that 
dynamical instabilities in the lengthening star streams 
in strong tidal fields are not likely to be present. For 
globular clusters, stream dyna mical instabilities are not 
seen in the n-body simulations ( Kiipper et al.|[2012l ). 

Cosmological n-body simulations from LCDM ini- 
tial conditions find that galactic halos will have 
approaching 10% of thei r mass in vast numbers 
of orbiting sub-halos (IDiemand. Kuhlen fe Madaul 
[20071: ISpringel et all [20081 iStadel et al l I2009D which 
will act on star streams to locally deflect and 



heat them (llbata et ail 
Sicgal-Gaskins & Valluri 




Johnston et all [2002; 
Carlbergj 120091) and 



induce visible gaps ( Yoon. Johnston fc Hogg) 120111 : 
ICarlberg||2012f) . The rate at which gaps are randomly 
created is approximately constant in time in a relatively 
slowly evolving, low redshift, galactic halo. Whereas 
the epicyclic pileups decrease with distance down the 
stream, sub-halo induced gaps increase in number and 
become better defined with age or distance down the 
stream. 

Here we take a dvantage of the SDSS DR8 data 
(|Aihara et al.l 120111 ) in the Pal 5 region to better ex- 



tract the density along the trailing stream. We use the 
adjoining sky regions to estimate the background and 
to estimate the errors in the density measurement. To 
characterize the gaps we develop a gap-finding filter. The 
resulting measurements of gap sizes and numbers are dis- 
cussed in relation to the ideas for epicyclic lumps and 
dark matter sub- halo induced gaps as a test of the LCDM 
dark matter sub-halo predictions. 

2. PAL 5 DATA ANALYSIS 



The SDSS Data Release 8 (|Aihara et al.l |20T1 has 
improved photometric uniformity over earlier releases 
along with a modest increase in sky coverage in the 
Pal 5 region. The stars are counted in 0.1° square 
pixels using the well - developed matc h ed fil ter technique 
ijRockosi et ail 120021 : iGrillmairl 120091 1201 ID . The pixel 
values are star counts filtered on the basis of their 
agreement with old, low metallicity stellar isochrones, 
weigh ted with the w Cen luminosity function (|de Marchil 
1999) using g, r and i photometry corrected for redden- 
ing. The binned, weighted star count data in the SDSS 
co-ordinate system are shown in Figure [T] 

2.1. The Extracted Stream 

The northern, trailin g, part of the stream is slig htly 
curved along its length (j Grillmair fc Dionatosl [20061 ) . as 
shown in Figure [TJ The stream centerline is defined first 
by rotating the stream to a horizontal x axis then remov- 
ing the residual variation with a spline function along 
the stream. We define the centerline of the stream with 
an eye fit, then iterate the positions but find that the 
improvement in the total integrated luminosity of the 
stream is not significant for a reasonable choice of the 
stream centerline. For each pixel we calculate xy co- 
ordinates relative to the centerline with the y (vertical) 
component set to zero. The pixels are maintained as sta- 
tistically independent measurement points for the anal- 
ysis. However, for the purpose of illustration, the data 
are interpolated back onto a square 0.1° grid, a 4° Gaus- 
sian smoothed image subtracted, and the difference is 
then smoothed with a 0.2° circular Gaussian. The two 
clusters are masked out. The image shown in Figure [2] 
places the fitted centerline at y=0 pixels, which leads to 
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Fig. 3. — The density profile transverse to the stream averaged 
over x=[40,272]. 

the cluster being position ed at fx,y1 = [4 2,-5] . At the dis- 
tance of Pal 5, 23.2 kpc (|Harrisl 11999 ) ■ each 0.1° pixel 
subtends approximately 0.0405 kpc. 

The stream is clearly detected over the entire length 
of the image, 42 to 274 pixels, an angular distance of 
23.2° from the cluster. The projected total length of the 
northern star stream is 9.4 kpc. The background rises 
slightly towards the end of the stream. We note that 
there is a negative density dip around x=215 pixels, with 
no clear source in the star counts or dust map. 

The bilinear interpolation and filtering used to create 
Figure [2] introduces correlations between the pixels, so 
all analysis is done using using the non-integer xy co- 
ordinates of the original pixels relative to the centerline. 
Figure [3] shows the density profile summed along the 
stream in 1 pixel bins transverse to the centerline. The 
FWHM of the stre am is approximately 3 pix els, equal 
to the 0.12 kpc that lOdenkirchen et all (|2003f ) estimate, 
but here over a nearly four times longer extent. We will 
normally do the analysis over a strip of total width of 
±2 pixels from the centerline, a full width of 4 pixels, 
capturing more than 95% of the stream density, which 
we use as the stream density profile in all of the analysis. 
The background is estimated from points on either side 
of the stream no closer than 5 pixels to the centerline, 
generally using two strips of five times the stream width, 
or 10 pixels. Using adjacent regions of 3 to 7 times the 
stream width changes the results about 10% or less. 

The two-dimensional map of Figure [5] shows that the 
stream has small deviations from t he centerline along its 
length , as was previously noted in iGrillmair fc Dionatod 
(2006), so some of the width is due to centerline wan- 
dering, not true width. We see no clear evidence for 
a systematic variation of stream width along its length. 
The full 2D map is potentially an extremely powerful 
tool in the analysis of the physical cause of variations in 
the stream density as the image signal to noise improves 
with better data. 

2.2. The Stream Density Profile and Its Errors 



Fig. 4. — The density along the stream (solid line) and in the 
surrounding background region (dotted line) in 0.5°, or 5 pixel, 
bins. The background region is 4 times the width of the stream to 
reduce noise. The cluster is at x=42. 
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Fig. 5. — The background subtracted density along the stream 
(solid line) and the error estimate (dashed line) in 0.5°, or 5 pixel, 
bins. 

The density along the stream and in the two adjacent 
regions used to estimate the background, binned over 5 
pixels along the stream to suppress noise, is shown in 
Figure [4] There is a slight rise in the background near 
the visible end of the stream, but it shows no small scale 
systematic deviations from the stream which would in- 
duce systematic errors that could be confused with gaps. 
The stream itself is detected at high significance along its 
entire length. The subtracted density profile along with 
the error estimate, both in the 5 pixels bins, is shown in 
Figure 

2.2.1. Random Errors 
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Constructing a stream density error estimate is crucial 
to the statistical tests for the significance of gaps. Our 
fundamental assumption is that the background level at 
the location of the stream and the stream error properties 
can be derived from the averages and standard deviations 
of the pixels in the regions to either side of the stream. 
The background itself is simply the average of the two 
sides. 

The standard deviation in a background pixel parallel 
to stream location i is calculated from pixels spread over 
a total width that is N times the width of the stream. 
Statistically, the background will have a standard de- 
viation of ai/y/N, where 0$ is the standard deviation 
in a region the width of the density stream. Accord- 
ingly, we take the error in the measurement of the derived 
stream density to be yl + N times the value measured 
in the background region, where we have made use of 
the fact that the stream is a relatively small, ~ 20%, 
over-density above the background, which in quadrature 
would change the errors about 10%. The exact value of 
<7j varies from pixel to pixel due to noise, but in practice 
it is convenient that there is no systematic trend along 
the stream, as shown in Figure [5j 

In detail, the stream density error estimate is con- 
structed using the background region a minimum of 5 
pixels from the centerline of the stream and generally 
five times the width of the stream on either side of it, all 
normalized to give the value appropriate for the stream 
itself. We have performed the analysis with varying back- 
ground regions from 3 to 7 times the width of the stream 
and see no differences above those expected from the es- 
timated noise level. The individual pixel density values 
in the background region, with a quadratic polynomial 
subtracted to remove the local mean, provides a basis 
for constructing gap-free bootstrap realizations of the 
stream density. That is, we construct random gap-less 
streams with a noise level the same as the real stream, 
using the background density profile along the stream, 
with the mean subtracted, scaled by \f\ + N and then 
create bootstrap samples using random drawing with re- 
placement. 

The background subtracted density in the 0.1° bins 
that we use for analysis is shown in Figure [6j along with 
a 0.2° Gaussian filtered version of the density to give 
better signal-to-noise for the purpose of the figure. The 
density profile in the first 6°, from x=40 to 100 recove rs 
the features first reported in IQdenkirchen et al.l (|2003f ). 

2.2.2. Extinction Systematic Errors 

Figure [7j shows the E(B-V) reddening values of 
ISchlegel et al.l ([1998D measured along the location of the 
centerline of the stream and the adjacent background re- 
gions. The reddening along the stream is typically slowly 
varying with a mean near 0.05 mag, requiring brightness 
corrections of 0.15 mag or less. The extinction in the re- 
gion we use for background is generally within 0.01 mag 
of the value along the stream. 

If the star counts with magnitude, m, were as steep as 
the Euclidean dlog (N)/dm oc 0.6m, an uncorrected uni- 
form distribution would have extinction induced density 
variations of ^23% at most, and more typically about 
7%. Because the extinction correction has been made 
and using the background to stream extinction differ- 
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Fig. 6. — The background subtracted density along the stream 
in 0.1° pixels (thin line) and filtered with a 0.2° width Gaussian 
(thick line). The Pal 5 cluster is at x=42. 
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Fig. 7. — The reddening along the stream (solid line) and in the 
background region (dashed line) in 0.1° pixels. 

ences of 0.01 mag or less as an estimate of the error the 
residual variations in number will be less than one per- 
cent, about a factor of ten less than the statistical errors. 
The reddening and extinction corrections introduce weak 
correlations between the pixels, since the originating data 
( Schl egel et al.lll998D are provided in 0.32° pixels, more 
than three times larger than the star count pixels. We 
conclude that any systematic error due to dust variations 
is less than 1%, hence is more than a factor of ten below 
the random errors. 

2.2.3. Random Error Consistency Test 

A simple x 2 test shows that the Pal 5 stream has den- 
sity variations far above the noise. The \ 2 P er degree of 
freedom for the data shown in Figure [6] gives \ 2 l v = 2-13 
for 224 degrees of freedom, which indicates very high sig- 
nificance variations, that is, less than a part in 10 6 that 
the density variations are due to noise . 

Along the stream we might expect the density of 
nearby points to be correlated. The auto-correlation 
of the mean subtracted density profile, 5d, as a func- 
tion of lag, cr 2 (Ax) = J S(x)8(x + Ax) dx, gives a cross- 
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Fig. 8. — The density profile as filtered with the wi (dotted thin 
line) and W2 (solid line) filters for gaps of width 5 pixels. Both 
curves have the same normalization, the standard deviation of the 
wi filter. 

correlation of adjacent 0.1° pixels of 0.2, which is not 
significant. Furthermore, the correlation analysis, even 
when restricted to the first 6°, reveals only statistically 
weak evidence for periodic structure, which would be ex- 
pected if the tidal cycloids dominated the density struc- 
ture along the stream. 

The auto-correlation function away from the initial 
peak provides an independent check on the size of the 
errors in the density distribution. The values of the cor- 
relation function away from the peak near zero lag are 
products of independent random variables whose popu- 
lation distribution has a standard deviation that is equal 
to the total variance of the data times the square root of 
the number of elements in the data. For the data range 
[50,274] the auto-correlation distribution gives the vari- 
ance per element to be 28.6. The measured x 2 /v = 2.13 
and the point-by-point estimate of the error over this 
range is 14.16. Therefore x 2 / v estimates the variance 
per element to be 30.1. The two independently deter- 
mined values agree within 5%, which requires that the 
independently estimated Ui values used in the \ 2 calcu- 
lation be correct. The agreement between two different 
methods to calculate the error and total variance gives us 
confidence in our overall error analysis. The next prob- 
lem is to assign some fraction of the total variance to a 
systematic signal of stream gaps. 

3. GAPS IN THE STREAM 

Sub-halo induced gaps will occur at random locations 
along the stream with the most distant, oldest, part of 
the stream expected to have relatively more gaps. Most 
gaps should be close in size to the smallest gaps that can 
survive the various blurring motions of the orbit, since 
the low mass sub-halos Narrow gaps should be more nu- 
merous because of the steep rise of sub-halo numbers 
with decreasing mass. Since mass is conserved, gaps 
in the streams must have approximately compensating 
positive and negative densities relative to a larger scale 
mean, although this assumption depends on measuring 
the density to sufficiently large distance from the stream 
centerline to include all scattered stars. The assumption 
will work best for smaller gaps. 



3.1. Gap Filters 

A procedure to find gaps is to convolve a function 
that approximates the expected shape with the den- 
sity data and identify peaks in the convolution. We 
use the functions w\{x) = (x 6 — 1) exp (—1.2321a; 2 ) and 
W2(x) = (x s — 1) exp (— 0.559x 4 ), which rise at x=0 from 
-1 through zero near x=l and then asymptotically back 
to zero, with integrals of (near) zero over [-3,3]. These fil- 
ters have horns on either side of a trough, approximating 
the shape of the density gaps found in the simulations of 
ICarlbergl (|2012f ) with W2{x) having a much flatter floor 
and sharper horns. The whole procedure is a form of 
wavelet analysis. As illustrated in Figure [5] the W2 filter 
generally separates close peaks better, so we preferen- 
tially quote the results it gives. 

We filter the density field with the two gap filters scal- 
ing them from 1 pixel up to 50 pixels in size and then 
find the peaks of the filtered field, corresponding to gaps 
in the density, and tabulate their heights and locations. 
To assign a confidence level to a given peak we run 
the same filters on our gap-less bootstrap samples con- 
structed from the background regions. The gap filters 
have the important property that they are symmetric 
with zero mean, so the mean value nor a linear slope of 
the data within a filter makes no difference to the out- 
come. Although not very important for these data, we 
use a low order polynomial fit to the background density 
to subtract the slowly varying mean background which 
would otherwise excessive variance in the bootstrap sam- 
ples. We usually generate 100,001 bootstrap samples us- 
ing the perl rand function, with randbit reported as 48. 
The heights of the peaks in the bootstrap samples are 
sorted and the heights which encompass, say 99%, of the 
sample, defines the 99% confidence levels as a function 
of the width of our gap filters. 

3.2. Gap Statistics 

The whole range of gaps and their confidence levels 
is shown in Figure [9] for the W2 filter scaled from 0.1° 
to 5.0°. The plot shows all gaps above 67% confidence. 
We note that most peaks are fairly stable, but that the 
peaks in the x=80 to 140 region merge together a set of 
narrow peaks into one very wide peak. The origin of this 
behavior is visible in Figure [6] where the density in this 
region becomes sufficiently low that the "one big gap" 
interpretation is plausible. However the confidence level 
assigned to this large gap is sensitive to the details of the 
filter shape as shown in the two panels of Figure [TO] 

To illustrate the outcome of the filtering we schemat- 
ically reconstitute the stream from the gap analysis. 
We sort the gaps by confidence level and insert gaps 
in a d(x) = 1 distribution with an illustrative depth 
of 0.5CL(90%)/CL at the location of the highest sig- 
nificance gap at any location. A lower significance gap 
within a higher significance gap is ignored. The outcome 
using gaps of 1.8° (750 pc) and smaller is shown in Fig- 
ure [TO] The wider filters merges a set of lower confidence 
peaks around x=120 into one. Although the graph is 
somewhat schematic the reconstituted gaps can be asso- 
ciated with features in the original data and with depths 
proportional to their statistical confidence. 

To count gaps we use the highest significance gap at 
any location ignoring any lower significance gap(s) within 
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Fig. 9. — The gap-finder peaks along the stream, plotted relative 
to the 90% confidence level for the W2 filter. 
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Fig. 10. — The gaps reconstructed from the highest significance 
gap filter peaks for gaps 1.8° and smaller (solid) and all gap sizes 
(dashed). The top plot is for the w\ filter and the bottom for 
the u>2 filter. See text for details. The plotted depth of the most 
significant gap at any location is 0.5 CL(90%)/CL, where CL is 
the confidence level. Only the 4 and 5 highest confidence gaps for 
filters 1 and 2, respectively, survive at 99% confidence The light 
line is the measured density of Fig. [6] 
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Fig. 11. — The stream density in bins of 4.7°. 

its range. In the case illustrated we count 14(10) gaps 
with 90% confidence or greater for the Wi(w2) filter. The 
region of negative density around x=215 may be a prob- 
lem in the star count analysis so, if present, we remove 
it from the peak count to reduce the counts by one. The 
fairly strong gap at x=55 is likely a tidal feature, dis- 
cussed below, that we also remove if present. There- 
fore, after discounting, we find at 90% confidence 12 (8) 
peaks. At 99% confidence we find 4(5) peaks. Moving 
to a very high confidence level, we find 1(2) peaks (after 
discounting) at 99.99% confidence. We will prefer the u>2 
filter results as better matched to the shape of gaps in 
our simulations and more stable counts. Our rates will 
conservatively use the 99% confidence counts, and use a 
multiplicative factor of two error, that is a range from 10 
to 2, with 5 gaps as the preferred value. 

3.3. Mean Stream Density 

The mean density of the stream over 4.7° is shown in 
Figure fTTI which illustrates the declining density with dis- 
tance from the progenitor. The decline is not likely to be 
a result of an increase of the tidal field, since at roughly 
7 Gyr old the Pal 5 stream has made several complete or- 
bits. The mean density decline could be due to a chang- 
ing mass loss rate from the cluster, but, simulations gen- 
erally expect a fairly constant of even declining mass los s 
rate with time (jKiipper et al.ll2010l: iDehnen et al.||2004fl . 
A relation to gaps is that since the most distant part of 
the stream is the oldest that material that was in the 
gaps moves sufficiently far from the stream that it is not 
captured in our density measurement which assumes a 
constant width along the stream. Of course it is quite 
possible that the observed stream exhibits several of the 
effects combined. 

4. GAP FORMATION PROCESSES 
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TABLE 1 

Observed Stream Gap Statistics 
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Fig. 12. — The gap rate-width relationship with all quantities 
scaled to 30 kpc radial distance from the galactic center. The 
dashed line is a variance weighted fit to the data. The dotted line 
is a model in which the number of halos varies as M -18 normalized 
to have the same number of M > 10 9 Mq sub-halos found for the 
M~ 19 fit of simulations. The hashed region is an estimate of the 
theoretical uncertainties for the cold-stream calculation. 

4.1. Tidal Density Variations 

iKiipper et al.l (|2008T ) find that the distance from the 
cluster to the first and subsequent epicyclic pileups is 
471-0, (4il 2 — k 2 )/k 3 , where xl is the distance from the 
cluster center to the Lagrange point and Vt and K are re- 
spectively the rotational and epicyclic frequencies. The 
Lagrange point is calculated as x\ — GM/(4f2 2 — k 2 ), 
where M is the cluster mass. For a locally flat rota- 
tion curve k = \J~2Q in which case the distan c e betw een 
epicyclic pileups is 2\/2 ttxl- IDehnen et al.l (|2004D es- 
timate the tidal radius with rf = GMR 2 /v 2 where R 
is the galactocentric radius and v c is the galactic cir- 
cular velocity at the cluster location, deriving r t to be 
0.054 kpc. The Lagrange point distance for a flat ro- 
tation curve is x\ = rf/2. Consequently xl = 0.043 
kpc for Pal 5. The epicyclic pileups at the base of the 
cycloidal type motion are then spaced 0.38 kpc apart 
which is 9.4 of our 0.1° pixels. These calculations apply 



to a circular orbit, whereas Pal 5's orbit is fai r ly ecc en- 
tric, for instance as modeled in IKiipper et all ()2012l ) or 
iMastrobuono-Battisti et all ()2012l ). but the same general 
considerations apply. In an eccentric orbit the spacing of 
the tidal stream density peaks and valleys vary as the 
stream moves around in the orbit. 

There are two types of tidal features worth seeking. 
The easiest are the lumps expected at the cluster po- 
sition, 42 pixels, plus the cycloidal distance, 9.4 pixels, 
or 51 pixels. There certainly is a density excess which 
peaks around 53 pixels, see Figures [2] and [6] There is a 
deficiency, not excess, at 60 pixels, but a clear excess at 
70 pixels. There are also deviations from the centerline 
in (our) y direction, the most prominent deviation being 
at about 32 pixels from Pal 5, or 1.3 kpc. There may be 
a similar deviation in the same direction at half the dis- 
tance, 16 pixels. These features seem sufficiently regular 
and of approximately the right character to be under- 
stood as the expected epicyclic pileups, but the distances 
are about a factor of two larger than the simple circu- 
lar orbit theory predicts, which could reflect the complex 
orbital and t idal history of this unusual low concentra- 
tion cluster (Kup per et al.l l2012). We discount the gap 
at x=58 pixels as more likely to be a cycloidal orbit gap 
than a sub-halo induced gap. 

4.2. Gaps and Sub- Halos 

Assigning the gap at x=58 to the epicyclic motion we 
remove that gap from our counts. Using the 99% con- 
fidence count of gaps we find 5 gaps. We estimate the 
age of the oldest part of the stream to be 7.2 Gyr, based 
on a drift rate at a velocity equal to the cluster disper- 
sion of 1.1 kms -1 and a length of the analyzed region of 
200 pixels, or, 8.1 kpc. The average age of the stream is 
therefore 3.6 Gyr. Consequently the average rate of gap 
creation is 0.17 gaps kpc -1 Gyr -1 for the Pal 5 trailing 
stream and we recommend a conservativ e error of a fac- 
tor of 2 on this rate. In an earlier paper (|Carlberdl20l"2T) 
we used the age of the oldest part of a set of streams to 
derive the gap rate. Modifying this to the average rate 
increases the gap rates a factor of two. Revised gap rates 
for the available streams using the mean age are given in 
Table 1. 

Sub-halos in galactic halos create gaps in a stellar 
streams at a rate which depends on the smallest visible 
gap which de pends on t he wi dth, w, of the stream, as 
worked out in[Carlberg (|2012l ). The earlier calculations 
were oriented towards structures further out in the halo 
and were scaled to 100 kpc. For this paper we undertook 
additional cold stream simulations that concentrated on 
the 15-30 kpc range and doubled the sub-halo orbit sam- 
pling. Since the gap width has a significant mass and 
radius dependence the simple fits to the results have a 
somewhat steeper slope with stream width, 



n u {w,r) = 5.0 x 10 



_ 4 n(r) 



"o 



r 0.45 



-1.16 



kpc 1 Gyr 



for w in kpc, r^o the orbital radius normalized to 30 kpc, 
and the constant evaluated for a stream age of 3.5 Gyr. 
For Pal 5's galacto centric distanc e of 19 kpc n(r) = 22n$ 
from Figure 11 of iSpringel et al.l ( 20081 ). For widths, w 
gr eater than a few kpc, Equation[T]is similar to the result 
of iCarlbersl (|2012D . but the steeper slope predicts rela- 
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tively more gaps in thin streams at small galactocentric 
radius than the same stream at a larger radius. 

The relation of Equation [T] is developed from a com- 
pletely cold, zero-width stream. The width in the re- 
lation arises from the assumption that the smallest gap 
length that is visible is the width of stream, on the basis 
that the random motions of the stars in the stream are 
well described by epicycles which have an extent perpen- 
dicular to the steam nearly equal t o the extent along the 
stream (|Binnev fc Tremaind 120081 ) . The cold stream as- 
sumption runs into a problem for low mass halos where 
the width of the stream becomes significantly greater 
than the scale radius of the sub-halo, beyond which 
stream stars are not sufficiently perturbed to create a vis- 
ible gap. For the Pal 5 stream of 0.12 kpc width and set- 
ting w = aR s (M) we find M w = 1.1 x 10 6 (l/a) 1 / - 43 M . 
The mass is close to the minimum mass halos that pro- 
duce the gaps and will lead to fewer visible gaps in 
warm stream simulations. A second diminishing effect 
is that there will be a small range in angular momenta 
(|Evre fc B inncv 2011), which translates into mean guid- 
ing center radius in the star streams such that differ- 
ential rotation will be an additional blurring effect that 
increases with distance down the stream. Consequently 
the cold stream rates will be upper limits to more de- 
tailed simulation results. 

The star stream gap-rate vs width relation appropri- 
ate for galactocentric radii less than 30 kpc is shown in 
Figure [T2] along with the available data. The dashed line 
is a variance weighted straight line fit. The fitted line 
has a slope of —0.59 ± 0.35, which is about 1.6 standard 
deviations below the predicted value. Adopting mean 
ages rather than end-point ages for the streams doubles 
the estimated gap rates and helps bring the data into 
agreement with the prediction. 



5. DISCUSSION AND CONCLUSIONS 

We have presented improved data which extends the 
known length of the northern (trailing) Pal 5 tidal stream 
to nearly 23°. The primary focus of this paper is to intro- 
duce an objective and statistically quantifiable approach 
to counting gaps, in the stream. At 90% confidence the 
results vary some 40% depending on the gap filter shape, 
but at 95% confidence and above the counts agree within 
one. At 99% confidence we find 5 gaps, where we have 
discounted one gap as being in a region in which the data 
are compromised. Our statistical confidence in the exis- 
tence of density variations in the stream is >99.9999% 
on the basis of its \ 2 - Not surprisingly, identifying the 
location and width of gaps can only be done at lower sta- 
tistical level, where we have settled on a relatively conser- 
vative value of 99% confidence. Deeper images would be 
invaluable in improving the statistics. Or alternatively, 
velocities for all of the current stars (g— 21 mag) would 
be invaluable in improving the statistical signal. 

The density structures within the first kiloparsec of 
Pal 5 are likely dominated by the orbital dynamics of 
stars escaping the cluster. However, at larger distances 
the rising number of gaps is consistent with the expecta- 
tions of dark matter sub-halo induced gaps. The small 
radius, cold stream, gap rate in terms of the minimum 



mass sub-halo M 8 = M/10 8 M Q is, 

Tl u (M 8 ,r) = 3.6 x lO-^rO^-^kpc- 1 Gyr" 1 , 

no 

(2) 

at 3.5 Gyr. The total number of sub-halos above the min- 
imum mass is 260Mg~ ' 9 . User the mid point and lower 
end of the stream averaged Pal 5 gap creation rate at 
r = 0.63 gives M$ = 0.15 - 1.0 x 10" 3 and implies a total 
population of some 1.3 — 7.2 x 10 5 sub-halos (over the 
433 kpc normalizing volume) or 900-5000 inside 30 kpc, 
far above the numbers of visibly populated dark mat- 
ter sub-halos in our galaxy. The inferred numbers are 
very sensitive to the rate estimates, varying as the 3.24 
power. An important limitation of our statistical count- 
ing of sub-halos is the use of completely cold stream in 
the simulations which give the most possible gaps at any 
sub-halo mass. The sub-halo masses we derive, about 
10 6 M Q are at the point where we expect that the stream 
width will reduce the number of persistent gaps. 

To gain an approximate understanding of the num- 
bers we take all the halos to have the same mass so that 
IZij ~ rnrbv, where n is their volume density, v is the 
typical encounter velocity, and b the largest distance of 
the typical (low mass) sub-halo that induces a gap. The 
value of b will be comparable to the stream width, £, so 
we set b = 1/2. For Pal 5 Ku ~ 0.17 kpc" 1 Gyr" 1 , 
and v ~ 220kms~\ We find n ~ 0.0041 kpc~ 3 . If 
these are the mean numbers inside 30 kpc, then the to- 
tal number inside this volume is about 1030 correcting 
the numbers for higher density at smaller radius with 
the Einasto mass function. Assuming a total mass inside 
30 kpc of 3 x 10 11 M Q and 7% of the mass in sub-halos 
gives a mean sub-halo mass of 2 x 1O 6 M That is, the 
huge numbers of halos required to create the gaps is a 
simple consequence of the fact that the interaction dis- 
tance for the smallest mass sub-halos that create gaps 
is comparable to the quite narrow stream width. The 
iSpringel et alJ (|2008h simulations find that sub-halos at 
these numbers and masses have a circular velocity near 
1 kms -1 , not surprisingly about equal to the velocity 
dispersion in the stream. It should be noted that for 
the mass function found in the simulations the heaviest 
sub-halos dominate the total mass. 

The main result of this paper is to put the counting 
of the gaps in the Pal 5 stream on more solid ground. 
The resulting gaps counts are reasonably statistically se- 
cure and their numbers supports the conclusion that the 
character and number of gaps in stellar streams is simi- 
lar to what a LCDM cosmological structure model pre- 
dicts. In comparison to a cold stream the number of gaps 
found is about a factor of three lower than predicted, but 
given the substantial uncertainties in both gap finding 
and in the simplified theory of gap creation used here, 
this mainly suggests that both the data and the theory 
n eeds to be made more precise. 

iCarlb crg (2001) stated: 11 If the lowest velocity disper- 
sion streams are older than about 3 Gyr, then LCDM 
subhalos would be ruled out\ noting that: 11 A firm con- 
clusion will require more extensive orbit modeling in 
well matched simulations." I Carlberg (|2012f ) provided ex- 
tensive but idealized circular cold stream modeling of 
sub-halo interactions with very low velocity dispersion 
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streams, reaching the conclusion that streams are eroded 
via gaps created at a linear rate at time. Figures HU1 and 
ITTI show that the Pal 5 stream has a fair fraction, say 
50% or so, of its stream removed in gaps. This stream is 
being eroded away and is harder to find as a continuous 
structure at large distance from the cluster. The exact 
fate of low velocity dispersion streams will be very de- 
pendent on their orbit around the galaxy. If they avoid 
dense regions such as the disk and especially the bulge, 
the stars pushed out of the gap remain near to the stream 
in a sort of snowplow pileup. However, if the stream or- 
bits close to the bulge then small orbital differences can 
quickly lead to wide physical separations. 

Recent discussion of pro blems with the number of 
the most massive sub-h alos ( Bo vlan-Kolchin et al.ll20~12t 
IStrigari &: Wec hslcr 2012|) has highlighted the sensitivity 
of the predicted sub-halo numbers to the overall mass 
normalization of the galaxy. A smaller overall galactic 
mass would somewhat improve the agreement of theory 
and observation for both the low and high mass end of 
the spectrum, but it will take additional physical process 
to resolve the discrepancy, one possibility being a ccount- 
ing for large scale flows ()Bovv fc Dvorkinl I2012T) . The 
concentration of sa t ellites and streams to a polar plane 
(|Lvnden-Belll [l~976t iPawlowski eFall 120121 ) is an addi- 
tional long standing puzzle that may call for a dwarf 
galaxy formation scenario with add itional astrophysics 
(|Hartwick|[200llKroupa et all 120 101) . We conclude that 
our indirect sub-halo counting technique, which relies 
only on the gravitational field of the sub-halos to induce 
star stream gaps, finds that the gap statistics requires a 
very large number of sub-halos, comparable to the num- 
bers that LCDM cosmology predicts. 

Ultimately the use of stream gaps to constrain the 
properties of density sub-structures within dark halos, 
that is, sub-halos, depends on being able to reliably count 



gaps, age date the tidal stream, and then relate the rate 
of gap creation to the population of sub-halos. In this pa- 
per we have made some progress on the statistics of gaps 
and are able to rely on 99% confidence gap detections. 
An unfortunate situation is that the inferred numbers 
are steeply dependent, the 3.2 power, of the gap creation 
rate. However, better stream modeling is entirely pos- 
sible and the observational statistics will improve with 
deeper images and as velocity data becomes available. 
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